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Wc study the back reaction of a thermal field in a weak gravitational background depicting the far-field limit of a 
black hole enclosed in a box by the Close Time Path (CTP) effective action and the influence functional method. We 
derive the noise and dissipation kernels of this system in terms of quantities in quasi-cquilibrium, and formally prove 
the existence of a Fluctuation-Dissipation Relation (FDR) at all temperatures between the quantum fluctuations of 
the thermal radiance and the dissipation of the gravitational field. This dynamical self-consistent interplay between 
the quantum field and the classical spacetime is, we believe, the correct way to treat back-reaction problems. To 
emphasize this point we derive an Einstcin-Langcvin equation which describes the non-equilibrium dynamics of the 
gravitational perturbations under the influence of the thermal field. Wc show the connection between our method 
and the linear response theory (LRT), and indicate how the functional method can provide more accurate results 
than prior derivations of FDRs via LRT in the test-field, static conditions. This method is in principle useful for 
treating fully non-equilibrium cases such as back reaction in black hole collapse. 



I. INTRODUCTION 



In a recent essay |l| one of us outlined the program of black hole fluctuations and back reaction we are pursuing 
using stochastic semiclassical gravity theory based on the Schwinger-Keldysh effective action || and the Feynman- 
Vernon influence functional [Q methods. We mentioned prior works for static, quasi-static and dynamic black hole 
spacetimes and commented on how to improve on their shortcomings. In this paper, following the cues suggested in 
that essay, we discuss first the static case of a Schwarzschild black hole, with a further simplificaton of taking the 
far-field limit, and derive the fluctuation-dissipation relation (FDR) || under these conditions. In our view (following 
Sciama's hints) the FDR embodies the back reaction of Hawking radiance . In a recent work || we treated a 
relativistic thermal plasma in a weak gravitational field. Since the far field limit of a Schwarzschild metric is just the 
perturbed Minkowski spacetime, the results there are useful for our present problem. The more complicated case of 
near-horizon limit is also doable by calculating the fluctuations of the energy momentum tensor for the quantum field 
(with the help of e.g., the Page approximation |^]). Derivation of the dissipation and noise kernels for a static black 
hole in a cavity are currently under investigation pO| , pT[ . 



A. Fluctuation and Back reaction in Static Black Holes 



We recapitulate what was said before [y on the state of art for problems in this case. Back reaction in this context 
usually refers to seeking a consistent solution of the semiclassical Einstein equation for the geometry of a black hole 
in equilibrium with its Hawking radiation (enclosed in a box to ensure relative stability). Much effort in the last 
15 years has been devoted to finding a regularized energy-momentum for the back-reaction calculation. (Sec []l2 13 
for recent status and earlier references.) Some important early work on back reaction was carried out by Bardeen, 
Hajicek and Israel JT3], and York [19], and more recently by Massar, Parentani and Piran |l6| along similar lines. 
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Since the quantum field in such problems is assumed to be in a Hartle-Hawking state, concepts and techniques from 
thermal field theory are useful. Hartle and Hawking Jl7j , Gibbons and Perry Q used the periodicity condition of the 
Green function on the Euclidean section to give a simple derivation of the Hawking temperature for a Schwarzschild 
black hole. The most relevant work to our present problem is that by Mottola (l9), who showed that in some generalized 
Hartle-Hawking states a FDR exists between the expectation values of the commutator and anti-commutator of the 
energy-momentum tensor. This FDR has a form familiar in linear response theory p0[: 
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where N and D are the anticommutator and commutator functions of the energy-momentum tensor, respectively (D 
is the temporal Fourier transform of D). That is, 

N a bcd(x,x') = ({f ab (x), f cd (x')})p 

D abcd {x, x') = ([f ab {x), f cd (x')])p. (2) 

He also identifies the two-point function D as a dissipation kernel by relating it to the time rate of change of the energy 
density when the metric is slightly perturbed. Thus, Eq.(|l|) represents a bonafide FDR relating the fluctuations of a 
certain quantity (say, energy density) to the time rate of change of the very same quantity. 

However, this type of FDR has rather restricted significance as it is based on the assumption of a specific background 
spacetime (static in this case) and state (thermal) of the matter field(s). It is not adequate for the description of back 
reaction where the spacetime and the state of matter are determined in a self-consistent manner by their dynamics 
and mutual influence. We should therefore look for a FDR for a parametric family of metrics (belonging to a general 
class) and a more general state of the quantum matter (in particular, for Boulware and Unruh states). We expect 
the derivation of such a FDR will be more complicated than the simple case above where the Green functions are 
periodic in imaginary time throughout (not just as an initial condition), and where one can simply take the results of 
linear response theory in thermal equilibrium (for all times) almost verbatim. 

Even in this simple case, it is worthwhile to note that there is a small departure from standard linear response 
theory for quantum systems. This arises from the observation that the dissipation kernel in usual linear response 
analyses is given by a two-point commutator function of the underlying quantum field, which is independent of the 
quantum state for free field theory. In this case, we are still restricted to free fields in a curved background. However, 
since the dissipation now depends on a two-point function of the stress-tensor, it is a four-point function of the field, 
with appropriate derivatives and coincidence limits. This function is, in general, state-dependent. We have seen 
examples from related cosmological back-reaction problems Q where it is possible to explicitly relate the dissipation 
to particle creation in the field, which is definitely a state-dependent process. For the black-hole case, this would 
imply a quantum-state-dependent damping of semiclassical perturbations. To obtain a causal FDR for states more 
general than the Hartle-Hawking state, one needs to use the in-in (or Schwinger-Keldysh) formalism applied to a class 
of quasistatic metrics (generalization of York [|l5j) and calculate the fluctuations of the energy mometum tensor for 
the noise kernel. In our problem such a calculation with back reaction is carried out in full detail, albeit only for a 
weak gravitational field here which depicts the far-field limit of a Schwarzschild black hole spacetimes. We wish to 
address the thermal field aspects of the problem, while saving the geometric aspects in the near-horizon case for a 
later investigation. 

B. Thermal Fields in Linear Gravity 

The behavior of a relativistic quantum field at finite temperature in a weak gravitational field has been studied 
before by Gross, Perry and Yaffe ^T|, Rebhan and coworkers [^2), de Almeida, Brandt, Frenkel and Taylor [^3| for 
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scalar and abelian gauge fields. In these work, the thermal graviton polarization tensor and the effective action have 
been calculated and applied to the study of the stability of hot fiat/curved spaces and "dynamics" of cosmological 
perturbations. To describe screening effects and stability of thermal quantum gravity, one needs only the real part 
of the polarization tensor, but for damping effects, the imaginary part is essential. The gravitational polarization 
tensor obtained from the thermal graviton self-energy represents only a part (the thermal correction to the vacuum 
polarization) of the finite temperature quantum stress tensor. There is in general also contributions from particle 
creation (from vacuum fluctuations at zero and finite temperatures). These processes engender dissipation in the 
dynamics of the gravitational field and their fluctuations appear as noise in the thermal field. We aim at finding a 
relation, the FDR, between these two processes, which embodies the back reaction self-consistently. 

In this work we use open system concepts and functional methods a la Schwinger-Keldysh j|] and Feynman- Vernon 
||. By casting the effective action in the form of an influence functional we derive the noise and dissipation kernels 
explicitly and prove that they satisfy a Fluctuation-Dissipation Relation (FDR) || at all temperatures. We also derive 
a stochastic semiclassical equation for the non-equilibrium dynamics of the gravitational field under the influence of 
the thermal radiance. 

We adopt the Hartle-Hawking picture where the black hole is bathed eternally - actually in quasi-thermal equilib- 
rium - in the Hawking radiance it emits. It will be described here by a massless scalar quantum field at the Hawking 
temperature. As is well-known this quasi-equilibrium condition is possible only if the black hole is enclosed in a box of 
size slightly larger than the event horizon jlB] (or embedded in an anti-de Sitter space Jm)). In the asymptotic limit, 
the gravitational field is described by a linear perturbation from Minkowski spacetime. In equilibrium the thermal 
bath can be characterized by a relativistic fluid with a four- velocity (time-like normalized vector field) u M , and temper- 
ature in its own rest frame /3 . Taking into account the four- velocity of the fluid, a manifestly Lorentz-covariant 
approach to thermal field theory may be used |2f|. However, in order to simplify the involved tensorial structure we 
work in the co-moving coordinate system of the fluid where it M = (1, 0, 0, 0). 

By making conformal transformations on the field and the spacetime, our results may be easily generalized to the 
case of a conformally coupled quantum scalar field at finite temperature in a spatially flat Friedmann-Robertson- 
Walker universe |2{|. Indeed we have earlier used the functional method and the Brownian motion paradigm Q to 
study similar problems in semiclassical gravity p7| . We found that quantum noise arising from fluctuations in the 
particle creation would constitute a stochastic source (whose effect can overdominate the expectation value of the 
energy momentum tensor in the semiclassical Einstein equation p8| ) in a new form of Einstein-Langevin equation [^) . 
We also came to the understanding that back reaction of vacuum quantum field processes (such as particle creation) 
on the dynamics of the early universe near the Planck time is summarily a manifestation of a FDR in semiclassical 
gravity. 

In Sec. ||, we describe our model and the derivation of the thermal CTP effective action. We compare it with 
the influence action Q and identify the dissipation and the noise kernels representing the linear response of the 
gravitational field and the quantum fluctuations of the thermal radiance, respectively. In Sec. Ill we show that they 
obey a fluctuation-dissipation relation at all temperatures. In Sec. IV we show how to derive a stochastic semiclassical 
equation for the gravitational perturbations from this effective action, which depicts the nonequilibrium dynamics of 
the gravitational field in a thermal radiation bath. 



II. CTP EFFECTIVE ACTION AT FINITE TEMPERATURE 
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A. The model 



In this section, we derive the CTP effective action for a thermal quantum field in a classical gravitational background. 
To describe the radiation we consider a free massless scalar field <p arbitrarily coupled to a gravitational field g^ v with 
classical action 

S m [<P,g^} = ~\J d n x ^[g^d^d„4> + t(n)R0 2 ] , (3) 

where R(x) is the scalar curvature and the arbitrary parameter £(n) defines the type of coupling between the scalar 
field and the gravitational field. If £(n) = (n — 2)/[4(n — 1)], where n is the spacetime dimensions, the field is said 
to be conformally coupled; if £(n) = the quantum field is said to be minimally coupled. In the weak field limit we 
consider a small perturbation from flat spacetime 7]^ v 

9n V {x) = T]^ + h^(x), (4) 

with signature (—,+,•••,+) for the Minkowski metric. Using this metric and neglecting the surface terms that appear 
in an integration by parts, the action for the scalar field may be written perturbatively as 

s m [(f>, VI = \ja n x<p[n + v« + v<v + ■■■]& (5) 

where the first and second order perturbative operators and are given by 
W = -\[d^{x)]d v + h^{x)d li d v +^(n)R {1 \x)Y 

= |[^^(x)]^ + ^(x)^^-e(n)(i? (2) (x) + ^(x) J R( 1 )(x))|. (6) 

In the above expressions, is the fc-order term in the pertubation h^ v {x) of the scalar curvature and the definitions 
hpv and h^ u denote a linear and a quadratic combination of the perturbation, respectively, 

hy, v = h*h au - ^hh^ + ^h 2 !]^ - ^h a[i h al3 -q^. (7) 
For the gravitational field we take the following action 



afL n - 4 
+ 4(n-4) 



J d n x ^—g 3R^ a0 (x)R^(x) - (l - 360(£(n) - l) 2 ) R(x)R(x) 



(8) 



The first term is the classical Einstein-Hilbert action and the second divergent term in four dimensions is the coun- 
terterm used in order to renormalize the effective action. In this action i 2 P — 16itG, a — (28807r 2 ) _1 and p, is an 
arbitrary mass scale. It is noteworthy that the counterterms are independent of the temperature because the thermal 
contribution to the effective action is finite and does not include additional divergencies. 

B. CTP effective action 

The CTP effective action at finite temperature for a free quantum scalar field in a gravitational background is given 

by 
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^ctpKA = Sf v [h+,]- Sf*[h- V ]- -Tr{h^[^]}, (9) 

where a, b = ± denote the forward and backward time path and G b [h^] is the complete 2x2 matrix propagator 

with thermal boundary conditions for the differential operator □ + + H . Although the actual form of & ab 

cannot be explicitly given, it is easy to obtain a perturbative expansion in terms of V Qb , the fc-order matrix version 
of the complete differential operator defined by V±± = ±V^ and V±^ = 0, and G b , the thermal matrix propagator 
for a massless scalar field in flat spacetime ||. To second order G b reads, 

U ab - U ab~ U ac V cd U db~ U ac V cd U db + U ac V cd U de V ef U fb H \ W ) 

Expanding the logarithm and dropping one term independent of the perturbation fv^ v (x), the CTP effective action 
may be perturbatively written as, 

r^pf/l^j = Sg^ [Jl^] — Sg W ^ 

+ l -Tr[vl l) Gl + V™Gt_ + vi 2) G p ++ - V {i) G f3 __] 



- l l Tr{vi 1) Gl + vi 1) Gl + + V PGt_vPGt_ - 2vf ) G^ 1) G^ + ]. (11) 



In computing the traces, some terms containing divergencies are canceled using countcrtcrms introduced in the 
classical gravitational action after dimensional regularization. In general, the non-local pieces are of the form 
Tr[Va G^V^ Gf s ]. In terms of the Fourier transformed thermal propagators G^ b (k) these traces can be writ- 
ten as, 

Tr[VpGi n V b W G? s \ = Jftft'^W^^') J * (2 J w e ifc -^gg tw (fc + g)G^(g)T^^( g ,fc), (12) 

where the tensor T fiV ' a0 (q, k) is defined in [0] after an expansion in terms of a basis of 14 tensors |p2]] . In particular, 
the last trace of ( |TT| ) may be split in two different kernels N A " /: "^(x — x') and D MI/,a/3 (x — x'), 

l -Tr[vl l) G + _v[ l) G _+] = - J d^x^x 1 h+ v {x)h- p {x')[T)^ a0 {x-x')+iW v ' a0 {x-x')]. (13) 

One can express the Fourier transforms of these kernels, respectively, as 



N (k) — J {9(k° + q o )0(-q°)+6(-k°-q o )9(q°)+n /3 (\q°\)+np(\k° + q°\) 

+2np{\q°\)np{\k° + q°\)} S(q 2 )6{(k + qf]T^(q, k), (14) 

D^ Q/3 (fc) = -in 2 J -^jL {9(k° + q°)8(- q °)-9(-k°- q °)8( q °)+ S g(k° + q °)np(\ q °\) 

- S g(q°)n (\k° + q°\)} S(q 2 )S[(k + q)*]T^(q, k). (15) 

Using the property T^' a/3 (<j, k) = T^ al3 {- q , -k), it is easy to see that W v ^{x-x') is symmetric and D A " y ' Q/3 \x-x') 
antisymmetric in their arguments; that is, W' af3 (x) = W> al3 (-x) and D^' Q/3 (a;) = -D WQ/J (-x). 

To properly identify the physical meanings of these kernels we have to write the renormalizcd CTP effective action 
at finite temperature ([□]) in an influence functional form Q]. N, the imaginary part of the CTP effective action can 
be identified with the noise kernel and D, the antisymmetric piece of the real part, with the dissipation kernel. In 
Sec. Ill we will see that these kernels thus identified indeed satisfy a thermal FDR. 

If we denote the difference and the sum of the perturbations h^ v defined along each branch C± of the complex time 
path of integration C by \h^ v ] = h~t v — h~ and {h^} = ht v + h^, respectively, the influence functional form of the 
thermal CTP effective action may be written to second order in as, 
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2£j 
1 
2 



-i / d 4 x dV [VlWD^fx-!')!^}^') 
+ J / AA' [h^](x)N^(x-x')[h a0 }(x'). 



(16) 



The first line is the Einstein-Hilbert action to second order in the perturbation h^ v (x). L^' (x) is a symmetric 
kernel (i.e. L^' a ^(a:) = L^ Q/3 (— x)) and its Fourier transform is given by 



L^TV) = \ [-k 2 Tr afl (q, k) + 2fc 2 Tf^(g, fc) + T^(q, k) - 2T^(q, k) 



(17) 



The fourteen elements of the tensor basis (q, k) [i = 1, • • • , 14) are defined in [^]. In the second line has 

the form of a perfect fluid stress-energy tensor 



rp\iV " 

" 30/? 4 



~(?f !/ + u / V) 
3 



(18) 



where is the four-velocity of the plasma and the factor is the familiar thermal energy density for massless 
scalar particles at temperature /3" 1 . In the third line, the Fourier transform of the symmetric kernel H^^^x) can 
be expressed as 



ah 



fc 4 n, \k 



H"(fc) = -^-^ln^Q^(fc) + ^Q— (fc) 



4 12 [i 



+^4 {-Tr a/3 («, k) 2T^(u, fc) + Tr Q V fc) + 2Tf'<*V k)} 



+ 96^ {^^'"V fc ) " 2fc 2 Tr Q/ V, fc) - Tf^^, fc) + 2T^ a/3 (u, fc)} 



+7T 



d 4 g 
(2^ 



<y(g 2 )^(|g°|)7' 



(fc + «z) 2 



+ ^[(fc + g) 2 ]n /3 (|fc° + g°|)P 



T^' Q/3 ((7,fc), (19) 



where p is a simple redefinition of the renormalization parameter \i given by /i = /iexp(|| + ^ ln47r — ^7), and the 



tensors Q^ v ' al3 (k) and Q^'" 13 (fc) are defined, respectively, by 

cr^(fc) = ^| T r Q/3 (g,fc)-^Tr Q/3 (g,fc) + ^T^ /? ('z,fc) 



[1 - 360(£ - ^) 2 ] Tr a/3 (9, fc) + ^T$r"( ? , fc) - ^T™(g, fc) !> , 



fc 4 



fc 2 



(20) 



Q^ Q/3 (fc) = [1 + 576(£ - i) 2 - 60(C - ~)(1 - 36£')] {t^ V fc) + ^T^fa fc) - ^T^fa fc)} . (21) 

In the above and subsequent equations, we denote the coupling parameter in four dimensions £(4) by £ and conse- 
quently £' means d£(n)/dn evaluated at n = 4. H M ' (fc) is the complete contribution of a free massless quantum 
scalar field to the thermal graviton polarization tensor |^,^3| and it is responsible for the instabilities found in flat 
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spacetime at finite temperature |H|~p3"|]. []. Eq. ( |l9| ) reflects the fact that the kernel & MU ' a ^ (k) has thermal as well as 
non-thermal contributions. Note that it reduces to the first term in the zero temperature limit (/3 — ► oo) 

*T*{k) ^ -^{li^q^(k) + lq^(k)Y (22) 

and at high temperatures the leading term (J3~ 4 ) may be written as 

^■ a0 (k) * ^fH.WTrKif), (23) 

where we have introduced the dimensionless external momentum = k^ /\k\ = (r, k). The H,-(r) coefficients were 
first given in and generalized to the next-to-leading order (/3~ 2 ) in ^3|. (They are given with the MTW sign 
convention in ||) 

Finally, as defined above, W' a(3 (x) is the noise kernel representing the random fluctuations of the thermal radiance 
and D^ ,a ^(x) is the dissipation kernel, describing the dissipation of energy of the gravitational field. 

III. FLUCTUATION-DISSIPATION RELATION AND LINEAR RESPONSE THEORY 

A. Fluctuation-Dissipation Relation 

These two kernels found above are functionally related by a fluctuation-dissipation relation (FDR). This relation 
reflects the balance between the quantum fluctuations in the thermal radiance and the energy loss by the gravitational 
field. In Q we have shown explicitly how this relation appears at zero temperature and at high temperature. Here 
using the properties of the thermal propagators, we show that the FDR is formally satisfied for all temperatures. 

We begin by showing that the FDR naively appears at the level of propagators as a direct consequence of the KMS 
relation ^9|. Then, using a generalization of this KMS relation, we see how the FDR is also satisfied by our noise 
and dissipation kernels. 

To obtain the FDR at the level of propagators we need to introduce the Schwinger and the Hadamard propagators. 
These propagators are defined as the thermal average of the anticommutator G(x — x') = —i([<f>(x), <j)(x')])p and the 
commutator G^\x — x') = — i({(fi(x), <p{x')Y)Pi respectively. The first represents the linear response of a relativistic 
system to an external perturbation and the second the random fluctuations of the system itself [^9],^) . Since we can 
write the KMS condition satisfied by the propagators G + _ and G_ + in Fourier space as 

&_ + {k) = e f3k °& + _{k), (24) 

the Fourier transform of both the Schwinger G{k) and the Hadamard Gp (k) propagators can be expressed, for 
example, in terms of G+_(fc) alone 

G(k) = &_ + {k) - Gi_{k) = (>° - l) & + _(k), 
Gf{k) = &_ + {k) + Gl_(k) = (e"*" + l)G£_(fc). (25) 

The FDR satisfied by these propagators follows inmediately from the above equalities tl~ 



1 Note that the addition of the contribution of other kinds of matter fields to the effective action, even graviton contributions, 
does not change the tensor structure of these kernels and only the overall factors are different to leading order p2[ 
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Gf{k) = cofo(Pp)G(k). (26) 

Obviously, this relation can also be recovered if we write the explicit expressions for the Fourier transform of the 
propagators 

G{k) = -2m sg{k°)5{k 2 ), 
Gf{k) = -2m coth (B^Tjs(k 2 ). (27) 

To use this last approach in our case could be a very difficult task because one needs to compute the integrals for the 
noise and dissipation kernels explicitly. On the other hand, if we follow the first technique we only need to generalize 
the KMS condition of Eq. ( p4] ) to the product of two propagators. This generalization reads 

& + {k + q)& + _{q) = eP k °& + _(k + q )& + (q), (28) 
and can be used to deduce the following formal identity 

& + (k + q)G + _(q)+& + _(k + q)&_ + (q) = coth f^j [& + {k + q)& + _{q) - & + _(k + q)& + {q)]. (29) 

Finally, one only needs to write, from the trace of Eq. ([l2]) and the definitions ( |l4| ) and ([15]) , the noise and dissipation 
kernels in terms of the propagators G±^, respectively, as 

N^V) = -i/(^4 [<?i+(fc + Q)Gi-(Q) +& + _{k + q)& + {q)]T^{q,k), (30) 

D^' Q/3 (fc) = 1/(04 [& + (k + q )GUq)-GUk + q)G^q)]T^(q,k), (31) 
and use the formal equality ( p9|) to prove that they are related by the thermal identity 

N^"V) = zcoth(^D^ Q/3 (fc). (32) 
In coordinate space we have the analogous expression 

N^(i) = J d 4 x' K FD (x- x')D^^(x'), (33) 
where the fluctuation-dissipation kernel Kfd(x — x') is given by the integral 

K FD (x - x>) = ij ^ e ik <*-* ') coth . (34) 

The proof of this FDR at finite temperature is in some sense formal because we have assumed along the argument 
that the integrals are always well defined. Nevertheless, the exact results obtained for the zero and high temperature 
limits Q indicate that the noise and dissipation kernels are well defined distributions p0[ |. The asymptotic analysis 
has also been useful to determine the physical origin of the fluctuations. 

B. Linear Response Theory 

A fluctuation-dissipation relation is usually derived using linear response theory (LRT). We now show the connection 
between the LRT [E0| and the functional methods we have used here. In the spirit of LRT the gravitational field is 
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considered as a weak external source which imparts disturbances to the radiance whose response is studied to linear 
order. 

Let us first recall the main features of LRT. Consider a system described by the Hamiltonian operator H Q initially 
coupled linearly to an external driving agent, say A a . Since we are only interested in how the system responds to 
the external agent, and not the details of the agent, we will ignore the Hamiltonian for the external perturbation but 
write the complete operator Hamiltonian of the sytem as 

H = H Q + A a J a , (35) 

where J a is the current operator associated with the external agent. If the system is in thermal equilibrium before 
the external source is applied the first order response of the system to this external force is given by the thermal 
expectation value of the commutator of the current operator over its thermal average 

([J»(x) - (J»(x))p, J v {x') - (J v {x')) ])p. (36) 

In contrast, the intrinsic quantum fluctuations of the system are described by the thermal average of the anticom- 
mutator. In our case, the conserved current operator is given by the stress-energy tensor T lll '{x) as derived from 
the classical action. Our objective is to show that the response and fluctuation functions for the stress-energy tensor 
considered in the LRT are equivalent to the dissipation and noise kernels, respectively. 

First, we write the classical action for the matter field to linear order in the gravitational perturbations, 

1 
2 

with the stress-energy tensor given by 

T» v = P^' a ^d a cf)df}(t> + C(v' iVa -d' 1 d u )<t> 2 . (38) 

Note that T^ v (x) is conserved if the classical unperturbed equation of motion for <fi is satisfied and it reduces to 
the stress-energy tensor for a scalar field in flat spacetime if £ = 0. Alternatively, we can write the Hamiltonian 
formulation of our problem. If we introduce the conjugate momentum variable of the matter field to first order in the 
perturbation 

dC -1, dT^ 

— r ~ <P+T:h uu — 

dcj) 2 d</> 



SUM/u,] - 7T / d 4 x + VTH, (37) 



n = — ~ HA"— ( 39 ) 



the Hamiltonian can be written as 

n 2 + (v^) 2 - • (40) 



(T^(x)T^(x')) - (T^(x)) (T^(x')) = -2 / e ik < x ~^ / & + (k + q)&l_(q)T^ (q, k). (41) 



Note that to first order cb and n are interchangeable in the expression for the stress-energy tensor. 

Using the thermal version of the Wick theorem |u],[32]], one can write, after some algebra, the equilibrium thermal 
average of the two-point function for the stress-energy tensor operator at different spacetime points in terms of 
products of thermal propagators, 

d 4 k ik-(x-x') f d 4 q 
(2tt)4 J (2tt)4 

Finally, defining A^T MI/ (x) = T^ u (x) — (T /JI/ (a;)) / 3 and using the expressions for the noise and dissipation kernels given 
in (30) and ( pl| ) respectively, we obtain 

{{A p T^{x),^T^{x')})p = 8N^(i-!'), (42) 

([A^O^A^V)]^ = ^V^ aP {x-x'). (43) 

From these formal identities we conclude that the functional method gives a description of the lowest order dynamics 
of a near-equilibrium system equivalent to that given traditionally by the LRT. 
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IV. EINSTEIN-LANGEVIN EQUATION 



To reinforce the points made at the beginning about the dynamic nature of the back reaction of thermal radiance 
on the black hole spacetime even for the quasi-static case, we now derive from the thermal CTP effective action a 
dynamical equation governing the dissipative evolution of the gravitational field under the influence of the fluctuations 
of the thermal radiance. It is in the form of an Einstein-Langevin equation ||. 

We first introduce the influence functional Q] T = exp(iSjF) where the influence action SVf is related to the CTP 
effective action in the semiclassical limit by 0, 



T = expz Ue{T CTp [h%]}+ 1 - J d 4 x dV [h^](x)N^(x - x')[h af3 ](x')\ , 



(44) 



where Re{ } denotes taking the real part. Following Q we can interpret the real part of the influence functional as 
the characteristic functional of a non-dynamical stochastic variable j^ v {x), 

$([V]) = exp J d 4 x d 4 x' [h^](x)N^(x - x')[h aP ](x')^j . (45) 

This classical stochastic field represents probabilistically the quantum fluctuations of the matter field and is responsible 
for the dissipation of the gravitational field. By definition, the above characteristic functional is the functional Fourier 
transform of the probability distribution functional V\j^ v \ with respect to j 1 *", 

$([VD = J Vf v V[f y ]e^ dix[h ^ ]{x)3 ' lV{x) . (46) 

Using ( [451) one can easily see that the probability distribution functional is related to the noise kernel by the formal 
expression, 



exp (-if d 4 x d 4 x' j^(x)[W u ' a P(x - x' 


T l M3(x>)) 


JVj^ exp(-±Jd 4 x d 4 x> j^{x)[N^ af3 (x 


- x')]-^ ]aP {x')) 



W] = 7 " T- ( 47 ) 

JVj^ exp(-±fd*x d*x> j^(x)[W^ al3 (x-x')}-^ a0 (x')j 

For an arbitrary functional of the stochastic field £[j^), the average value with respect to the previous probability 
distribution functional is defined as the functional integral (£[j^])j = J T)\]^ v \ V\j^ v \E In terms of this 
stochastic average the influence functional can be written as T = (exp (i^cTp[^uiA))ji where ^ctpI^uu] ^ s the 
modified effective action 

^ctpKA = Re{T? CTP [h±„]} + J d 4 x [h^]{x)r{x). (48) 

Clearly, because of the quadratic nature of the characteristic functional ( f45| ) and its relation with the probability 
distribution functional (f46|), the field j^ix) is a zero mean Gaussian stochastic variable. This means that its two- 
point correlation function, which is given in terms of the noise kernel by 

(r{x)3 ap (x')) 3 = N^(x-x'), (49) 

completely characterizes the stochastic process. The Einstein-Langevin equation follows from taking the functional 
derivative of the stochastic effective action (|4g| ) with respect to [hfj, v ](x) and imposing [h^ u ](x) = 0. In our case, this 
leads to 



1 



-pr Id x L 



4 / t fiv 7 a0 



-P ' 



^(x-x')h a(3 (x') + ±T^+ Jd 4 x' (K^(x-x')-n^ aP {x-x'))h a p{x')+r{x) = 0. (50) 



To obtain a simpler and clearer expression we can rewrite this stochastic equation for the gravitational perturbation 
in the harmonic gauge W v v = 0, 
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Uh^{x) + £% + 2P paMj j dV (K^(x - x') - B^ a(3 (x - a;')) + 2^(^)1 - 0, (51) 

where we have used the definition for h^ix) written in (0) and the tensor P pa . a is given by 

Ppo,aP = ^ (VpaVcrP + llpPVaa - VpaVap) ■ (52) 

Note that this differential stochastic equation includes a non-local term responsible for the dissipation of the gravita- 
tional field and a noise source term which accounts for the fluctuations in the thermal radiance. They are connected 
by a FDR as described in the last section. Note also that this equation in combination with the correlation for the 
stochastic variable (^9|) determine the two-point correlation for the stochastic metric fluctuations (h^ (x)h a p(x'))j 
self-consistently. 

V. CONCLUSIONS 

In this paper we show how the functional methods can be used effectively to study the non-equilibrium dynamics 
of a weak classical gravitational field in a thermal quantum field. The Close Time Path (CTP) effective action and 
the influcence functional were used to derive the noise and dissipation kernels of this system. The back reaction of the 
thermal radiance on the gravitational field is embodied in a Fluctuation-Dissipation Relation (FDR) , which connects 
the fluctuations in the thermal radiation and the energy dissipation of the gravitational field. We prove formally the 
existence of such a relation for thermal fields at all temperatures. 

We also show the formal equivalence of this method with Linear Response Theory (LRT) for lowest order perturbance 
of a near-equilibrium system, and how the response functions such as the contribution of the quantum scalar field 
to the thermal graviton polarization tensor can be derived. An important quantity not usually obtained in LRT 
but of equal importance manifest in the CTP approach is the noise term arising from the quantum and statistical 
fluctuations in the thermal field. 

Finally, we emphasize that the back reaction is intrinsically a dynamic process which traditional LRT calculations 
cannot capture fully. We illustrate this point by deriving a Einstein-Langevin equation for the non-equilibrium 
dynamics of the gravitational field with back reaction from the thermal field. This method can be applied to quasi- 
dynamic H[53| or fully dynamic problems such as black hole collapse ]l^,|l6|,^4) . To complete the present problem of 
quasi-static black hole back reaction, we need to perform the same calculation for the full Schwarzschild spacctimc. 
Currently we are working on the fluctuations of the energy-momentum tensor near the black hole horizon and the 
derivation of the noise kernel. Results will be reported in future publications. 
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